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Abstract 

Measurements of Pitot pressure were 
made in the exit plane and plume of a low- 
density, nitrogen nozzle now. Two numerica 
computer codes were used to analyze the 
flow including one based on continuum theory 
using the explicit MacCormack method, and 
the other on kinetic theory using the m< Bthod 
of direct-simulation Monte Carlo (DS )• 
continuum analysis was carried to the nozzle 
exit plane and the results were compared to 
the measurements. The DSMC analysis was 
extended into the plume of the nozzle flow a 
the results were compared with measure 
ments at the exit plane and axial stations 12, 
24 ar,d 36 mm into the hear-field p »m Two 
exDerimental apparatus were used that dit 
fered in design and gave slightly different 
profiles of pressure measurements^ 

DSMC method compared well with the mea 
urements from each apparatus at all axial 
stations and provided a more accurate p 
diction of the flow than the continuum method, 
verifying the validity of DSMC for such calcu- 
iations. 


Introduction 

NASA has a continuing effort in developing 
small rockets that operate on electric power 
for both primary and auxiliary propulsion on 
satellites. These rockets are typ'ned by 
thrusts of order 1 N or less and also by 
viscous nozzle flows because of their small 
size and somewhat high operating temper- 


atures. An important consideration in i Pacing 
such thrusters on a satellite is the effec 
the plumes may have on satellite SUI t a ?® s ® 
systems. Given the viscous nature of the flow, 
the plumes tend to spread widely in vacuu 
witha considerable now in the backward di- 
rection Some major effects from the plum 
could include contamination from mass depo- 
sition and unsatisfactory disturbance torques 
or thermal loading from impingement on van 
ous satellite surfaces. 

To gain a better understanding °f thrus J® r * 
satellite interaction and design considerations 
in placing electric propulsion on satellites, 

• study of thruster plumes is currently in 
progress. Assessment of plume effects from 
small thrusters is being approached analyt 
icallv by modeling the nozzle flow and plume 
on "L mSecular level. Laboratory lea s are 
also beinq conducted in a vacuum facility to 
complement the analyses by P r ° v| d'"9 ®* p ® ’ 
imental measurements for verification of the 
numerical codes. Under study are low 
density, viscous nozzle flows and the P |u mes 
that these flows generate in a vacuum T 
flows are characteristic of low-thrust device 
such as arcjets and resistojets. Of particular 
interest is the prediction of P'^me expansion 
especially in the off-axis region, that may a 
feet surfaces in the proximity of a thruster. 

In a previous study, the flow of carbon 
dioxide in a nozzle was computed with two 
numerical techniques predicated on different 
flow regimes 1 . The first technique was based 
on continuum theory and numerically solve 
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the Navier-Stokes equations for compressible 
flow. The second technique was based on 
kinetic theory and used the method of direct 
simulation Monte Carlo (DSMC). Both nu- 
merical methods were applied to a reference 
problem defined to be a low-density, viscous 
gas flow in a convergent-divergent nozzle of 
conical shape. 

The work described in Reference 1 dem- 
onstrated that the numerically intensive DSMC 
technique could be readily applied to a low- 
density nozzle flow, where the flow varied 
from continuum (Kn~10- 3 ) at the throat to 
rarefied (Kn~10 -1 ) at the exit plane, where 
Kn is the Knudsen number based on local 
nozzle diameter. The solutions from the two 
techniques exhibited differences, mainly in the 
region near the wall at the nozzle lip In a 
subsequent study 3 , nitrogen was used in place 
of carbon dioxide to remove the ambiguity of 
specifying an effective specific heat ratio in 
the continuum method, and allowed for a 
more direct comparison between the analyt- 
ical schemes. 

To provide verification of the analyses, the 
reference problem was also run as an exper- 
iment in a large vacuum facility. The objective 
of the experiment was to measure pressure in 
the nozzle exit plane and near-field plume 
with Pitot tubes Similar work was done by 
Bailey and Price 3 for flow of C0 2 , and served 
as a guide for the work of this study. This 
paper presents pressure measurements taken 
in a nitrogen nozzle plume and the compara- 
ble pressures computed from the two numer- 
ical codes. 


Reference Problem 

The reference problem was defined to typ- 
ify the nozzle flow of an electric propulsion 
device while also matching the flow capacity 
of the vacuum pumps on the experimental fa- 
cility. A conical nozzle geometry was se- 
lected for both analytical and machining 
simplicity, but it is also typical of small 
thruster nozzles. The geometry for the prob- 
lem is illustrated in Figure 1, with details given 
in Table 1. The nozzle operating parameters 
are listed in Table 2 for two configurations 
corresponding to two different experimental 
apparatus and two operating conditions. 

The nominal Reynolds Number of 850, 
based on throat diameter and stagnation 
viscosity, is characteristic of a high- 


temperature, high-performance resistojet 4 . 
The nozzle was made sufficiently large to 
minimize probe size relative to the exit diam- 
eter and does not resemble actual thruster 
hardware. Nitrogen was selected as the test 
gas for analytical convenience and, similarity, 
does not particularly represent an actual 
working fluid of a spacecraft propulsion sys- 
tem 


Test Facility 

The experimental tests for this study were 
conducted in a space-simulation vacuum tank. 
The test apparatus was mounted in a 1 m di- 
ameter section about 1 m in length, and at- 
tached to a larger tank 4 9 m in diameter and 
19 m long. The facility pumping system con- 
sists of twenty oil-diffusion pumps in series 
with four blowers and four roughing pumps 
A detailed description of the test facility can 
be found in Reference 5. 

The pumping system was able to maintain 
a vacuum during the tests of order 10 " 2 Pa for 
a nozzle flowrate of 6.8 x 10- 5 kg/s. The vac- 
uum pressure was monitored with a hot- 
cathode ionization gage mounted on the test 
section and connected to a digital meter 


Experimental Apparatus 

The test apparatus consisted of a heat 
exchanger and nozzle that simulated a 
thruster and traversing tables that permitted 
surveying the plume with pressure probes. 
Two heat-exchanger assemblies of different 
design were used in the tests. The purpose 
of the second design was to reduce the tem- 
perature of the heating element to lengthen its 
operational life. 

In the original apparatus {configuration 1), 
the nitrogen was heated in a 3.2 mm diameter 
tube that was coiled around a 15.9 mm car- 
tridge heater. In the newer apparatus (con- 
figuration 2). the nitrogen was heated in an 
annular area comprised of a 12.7 mm diam- 
eter cartridge heater contained in a larger 
tube of 17.3 mm internal diameter. A sche- 
matic of this configuration is shown in Figure 
2. The heater for configuration 2 had a lower 
power density than configuration 1. The heat 
exchangers produced the same gas temper- 
ature, but the second design lowered the heat 
conduction along the heat-exchanger and 
nozzle walls and produced nozzle wall tem- 
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sure scans were taken only in the horizontal 
plane starting at the nozzle centerline (R = 0). 


Analytical Methods 

Two methods of analysis were applied to 
the problem. One was a numerical code 
based on continuum flow, and the other a 
code based on kinetic theory that considered 
molecular motion by the method of direct- 
simulation Monte Carlo. 


Continuum Method 

The continuum computation of the flow was 
made with the numerical code termed the 

«l!i C A°on „ Nozzle Analysis Program, version 2 
(VNAP2) 8 . The code solves the Navier-Stokes 
equations for a compressible fluid in time- 
dependent, non-conservative form by the ex- 
plicit, two-step MacCormack method which 
gives second-order accuracy in space. With 
this scheme the equations are marched in 
time from a specified initial flowfield to a 
steady-state solution. For this problem, given 
the relatively low Reynolds number, laminar 
flow was assumed. Also assumed was an 
adiabatic nozzle wall. 

The code produces a body-fitted, non- 
orthogonal grid in physical coordinates that, 
by transformation, is rectangular in computa- 
tional space. For the problem in this study, 

1 radial and 51 axial grid lines were speci- 
fied. Furthermore, the grid was clustered in 
the nozzle throat and near the nozzle wall to 
capture adequately the steep gradients of the 
flow variables in those regions. 

The computation was started in the up- 
stream, subsonic portion of the nozzle The 
inflow conditions included total pressure and 
temperature, and zero radial velocity. The 
subsonic inflow boundary employed the 
method of characteristics such that the axial 
velocity and static pressure and temperature 
were determined from interior, downstream 
points, and evolved with the solution from ini- 
tial values. 


plane and was not known a priori, but did 
match the ambient at. some finite distance into 
the plume. This condition was simulated by 
extrapolating the exit-plane pressure from in- 
terior points (as normally done in the super- 
sonic flow) in the subsonic region at each time 
step and along each grid line, then using the 
new value as the specified ambient pressure 
In this manner the exit-plane pressure was 
continually updated iteratively. This scheme 
computed relatively smooth profiles of prop- 
erties in the vicinity of the wall at the exit 
plane, but did impose an effective back pres- 
sure on the flow. In cases where the subsonic 
exit-plane pressure was explicitly fixed, the 
properties exhibited a discontinuity across the 
sonic line where the flow decelerates from 
supersonic to subsonic in the direction normal 
to the wall. 

The scheme was marched for 20,000 
timesteps from an initial condition corre- 
sponding to one-dimensional, isentropic flow 
The solution required about 0 8 h of CPU time 
on a Cray-YMP. Convergence of the scheme 
was not monitored as such but conservation 
of mass was checked and used as a guide 
At 20,000 timesteps the mass flow was fairly 
consistent throughout the nozzle. In the di- 
verging section it was within 1% of the value 
specified by the inlet conditions, and in the 
converging section, within 2%. 


Monte Carlo Method 

The Monte Carlo method used to simulate 
the nozzle flow was developed by Boyd’ after 
the method of Bird*. The code was highly op- 
timized and specifically structured for 
vectorization on a supercomputer, and per- 
mitted the practical use of a large number of 
simulated molecules and computational cells 
necessary to model an internal flow. 


Since the stream exhausted to an ambient 
pressure near zero, the outflow boundary 
condition for the subsonic portion of the flow 
in the exit plane was modified from the ori- 
ginal scheme which required the explicit 
specification of an ambient pressure. It was 
assumed that the pressure of the subsonic 
flow did not match the ambient at the exit 


Particle collisions were simulated in the 
code by the Variable Hard Sphere model of 
Bird*. The number of collision pairs sampled 
in each cell is based on a probability function 
dependent on particle relative velocity and 
collision cross-section. Rotational nonequi- 
librium was included in the scheme with a re- 
lationship developed by Boyd 9 for the 
probability of energy transfer between the 
translational and rotational modes and was 
fully integrated into the vectorized procedures 
of the code. 
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Calculation of Pitot Pressure 


To compare the numerical solutions with 
experimental data, Pitot pressures were first 
calculated from the numerical results usinq 
the computed state variables. The calculated 
Pitot pressure can be viewed as the pressure 
that would be measured if a Pitot tube were 
inserted into the stream computed by the nu- 
merical codes. Since the flow was both 
supersonic and rarefied, the calculation in- 
volved a two-step process: 

1 The pressure ratio across a normal shock 
was calculated with the Rayleigh Pitot tube 
equation (c.f. Ref 11, p 154) represented by 

P oy!P, = f(y, M/) ( 1 ) 

where P oy is the total pressure behind the 
shock, P, and M„ are the static pressure and 
Mach number ahead of the shock obtained 
from the numerical solutions, and v = 1.4 j s the 
ratio of specific heats. 


2. The ideal pressure P oy was corrected 
rarefaction effects by 


for 


i ne viscosity, p y , was derived from the nu- 
merical results by first calculating the tem- 
perature T y . This quantity was computed from 
the normal-shock relation for static temper- 
ature (c.f. Ref. 11 , p. ns) represented by 


TylL = f( y , M/) 


(4) 


w At re ,l l® va,ues of static temperature ahead 
o the shock, T Xt and the Mach number, M 
were obtained from the numerical solutions 
The viscosity was then found from the stand- 
ard temperature power-law relating the 
viscosity of gases (c.f. Ref 13 p 28) 


Py Pref (Ty/TretT 


(5) 


where n rr , is the value for viscosity of nitrogen 
at the reference temperature, T, * 


n = 0.75. 


and 


.. ^ Presentation of analytical results from 
the DSMC solutions (except those of Figure 7) 
utilize P^. No corrections were made to the 
ritot-pressure measurements. 


lo 9( P o ml Poy) = 089 - .120 log(fle p ) (2) 


Experimental and Analytical Results 


Re p < 5.6 

In equation (2), Re p is the probe Reynolc 
number and P om is the corrected total pressur 
behind the shock and is the pressure thj 
would actually be measured with a Pitot tub( 
Equation (2) is a fit of data obtained from 
report on rarefaction effects of Pitot tubes fo 
measuring pressure in rarefied, supersonic 
wind-tunnel flow 12 . In Reference 12, the facto 
^cmlPoy, which is always > 1, was plotted as 
function of Re p , with total temperature of the 
stream, T p , as parameter Equation (2) is a fi 
of the parametric curve for T a = 700 K. 


In Reference 12, the 
ber was defined as 


probe Reynolds num- 


Re p = PooU„D p /fi y ( 3 ) 

w /\® re f°° and U °° are ,he der >sity and velocity 
of the free stream preceding the shock, D„ is 
the probe diameter , and p y , the viscosity of the 
gas at the stream temperature, T y , behind the 
shock. In the probe Reynolds number, the 
actual diameter of the Pitot tube on the ex- 
perimental apparatus was used for D p . The 
values of and were obtained directly 

output of the codes, except for the 
DSMC results, where the mass density was 
first calculated from number density the ac- 
tual quantity computed by the code 


Radial scans of Pitot pressures for both 
experimental configurations, at an axial sta- 
tion 12 mm from the exit plane of the nozzle 
are shown in Figure 6. The pressures from 
configuration 2 (triangular symbol) are lower 
and the profile flatter, in the core flow near the 
centerline than those for configuration 1 (cir- 
cular symbol). The pressures from both con- 
figurations are nearly the same in the outer 
region of the plume. This figure is typical of 
the difference in pressure profiles between 
the two experimental configurations. 

A comparison of results from the numerical 
methods and the experimental data at the 
nozzle exit plane (Z=0) is given in Figure 7 for 
configuration 1. Plotted in the figure is Pitot 
pressure, both measured and calculated from 
each of the numerical methods, as a function 
of radial distance. The radial scale starts at 
the nozzle axis and extends to the nozzle wall 
at 16 mm. In this figure the results from the 
numerical methods are not corrected for 
rarefaction effects, so the numerical results 
are those for the ideal, total pressure, P^, be- 
hind a normal shock. Because of the fairly 
large region of subsonic flow in the continuum 
solution, the correction for rarefaction effects 
could not be applied to a considerable portion 
of the flow (a subsonic probe-recovery factor 
would have to be used). Furthermore, in the 
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core region near the centerline where the 
variation in pressure is greatest between the 
two analyses, Re p > 5.6, making the rarefaction 
correction nil. Thus for a fair portion of the 
flow, Figure 7 presents a valid comparison of 
the pressures. It can be seen from the figure 
that the DSMC computation gives the better 
correspondence to the experimental data. 
The continuum computation, while yielding a 
higher pressure in the core, yields a slightly 
lower pressure in the viscous region near the 
wall. 

In Figure 8, Pitot pressures computed from 
the DSMC results {all remaining results are 
corrected for rarefaction effects) and exper- 
imental data are given at the axial station 12 
mm from the exit plane for configuration 1. 
The profiles of Pitot pressure extend radially 
from the nozzle centerline to 24 mm. From the 
figure, the DSMC solution is seen to give good 
agreement to the experimental measure- 
ments. 

Similar plots of DSMC results and exper- 
imental data for configuration 1, at axial 
stations of 24 mm and 36 mm from the nozzle 
exit plane are given in Figures 9 and 10. In 
these figures, the radial scan is also from the 
nozzle centerline to 24 mm. Again, the corre- 
spondence is quite good. 

Pitot pressures computed from the DSMC 
results and measured values for configuration 
2 are presented in Figures 11 through 14. The 
axial stations are, as before, at the exit plane 
(Z = 0) and 12, 24, and 36 mm downstream into 
the plume. In general, the DSMC solution 
reasonably predicts the measured pressure 
profiles. In Figures 12 and 13, at the 12 and 
24 mm axial stations, the centerline pressure 
from the DSMC solution is nearly identical to 
the measured value. The pressures from the 
DSMC solution are slightly lower at all radial 
positions for the 36 mm axial station, Figure 
14. 

Finally, Figure 15 shows a comparison of 
flow angles from the DSMC solution and 
measured values from configuration 2 for the 
axial station 12 mm from the exit plane. There 
is fair agreement between the computed and 
measured values of flow angle, with the larg- 
est descrepancy occuring at a radial distance 
of about 10 mm. The experimental measure- 
ments exhibited some scatter in the region 
from 8 mm to 12 mm in radial distance. 


Discussion of Results 

In general the correspondence between 
the DSMC results and experimental data is 
better for configuration 1, but the DSMC sol- 
utions do predict the difference in pressure 
profiles between the two experimental appa- 
ratus. As previously mentioned, the principal 
difference in experimental configurations was 
the design of the heat exchangers. The ear- 
lier design (configuration 1) produced nozzle 
wall temperatures about 40°K higher because 
of higher heat conduction along the tube of the 
heat exchanger to the nozzle wall. The cooler 
wall of configuration 2 would tend to flatten the 
pressure profile as the density and pressure 
of the flow near the wall would be higher rel- 
ative to the centerline region. Also, for the 
same flowrate, the total pressure of 6210 Pa 
for the nozzle of configuration 2 was slightly 
lower than the 6400 Pa of configuration 1 
which would give lower stream pressures 
overall for configuration 2. 

The continuum code was not run for the 
conditions of configuration 2, so properties of 
the inflow surface were not provided the 
DSMC model for this case. Rather, the inflow 
number density for the DSMC simulation was 
scaled from the perfect gas relation, where 
the stream static temperature was the same 
as the first case. Nozzle surface temperatures 
were scaled from experimental measure- 
ments to simulate the wall boundary condition 
for configuration 2. Although the DSMC sim- 
ulation for configuration 2 needs some refine- 
ment, the computed results, nonetheless, 
agree reasonably well with the experimental 
data. 

One factor that would affect the compar- 
ison between computed and experimental re- 
sults is the actual ambient pressure of the test 
facility In particular, the facility pressure 
during testing was about 10- 2 Pa whereas an 
absolute vacuum was essentially assumed for 
the DSMC simulation. The ambient pressure 
of the test facility would restrict plume expan- 
sion and have an effect on both pressure and 
flow angle, especially in the far-field regions 
of the plume. This may be the reason that the 
pressures from the DSMC simulation were 
slightly less than measured values at an axial 
distance of 36 mm as given in Figure 14. 
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Summary and Conclusion 

A nozzle flow typical of small thrusters 
used on satellites was simulated by an ex- 
perimental apparatus and Pitot pressures 
measured in the exit plane of the nozzle and 
at three axial stations downstream of the exit 
plane. The measurements were made on two 
experimental apparatus differing in heat- 
exchanger design and corresponding to two 
slightly different operating conditions. The 
measurements were made to specifically val- 
idate computed results from a kinetic theory 
model of the flow using the DSMC technique. 
The measurements were also used to check 
results from a continuum model, although the 
continuum model was used primarily to spec- 
ify an inflow surface and aid in defining a grid 
for the DSMC model. Comparisons were 
made between the measured pressures and 
the equivalent pressure derived from the nu- 
merical results, including a correction for 
rarefaction effects. 

Design differences in the laboratory appa- 
ratus resulted in differences in the measured 
profiles of Pitot pressure that were corrob- 
orated by the DSMC method when the appro- 
priate pressure and wall temperatures were 
used in the simulation. Agreement between 
the experimental measurements and the 
DSMC results was quite good which verifies 
the validity of the method for computing low- 
density nozzle and plume flows. The com- 
parisons emphasize the importance of having 
experimental data available from a laboratory 
apparatus for validating a numerical simu- 
lation. 
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Tabic 1. Nozzle Geometry 


Throat diameter, Dt 

3.18 

mm 

Exit 

diameter, D e 

31.8 

mm 

Inlet 

diameter, Dj 

22.1 

mm 

Wall 

thickness, t w 

1.65 

mm 

Area 

ratio, (D e /D t ) 2 

100 


Inlet 

half-angle, ©i 

45° 


Exit 

half-angle, 0 e 

20° 



Table 2. Nozzle Flow Conditions for the 

Two Experimental Configurations 


Configuration 

1 

2 

Total pressure, P 0 

6400 Pa 

6210 Pa 

Total temperature, T 0 

699 K 

700 K 

Mass flow, rh 

6.8xl0' 5 kg/s 

6.8xl0' 5 kg/s 

Reynolds number, Re* 

850 

850 

Wall temperature, T w i 

551 K 

507 K 

Wall temperature, T w 2 

539 K 

498 K 


* Re = 4m/D t p 0 , where |i 0 is the gas viscosity at T c 



Figure 1. Nozzle Diagram 
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Figure 6. Measured Pitot Pressure Profiles for Configuration 1 and 2 

(Axial Location, Z = 12 mm) 
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Figure 7. Comparison of Pitot Pressure Profiles for Configuration 1 

(Axial Location, Z = 0 mm) 
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Figure 8. Measured and DSMC Pitot Pressure Profiles for Configuration 1 

(Axial Location, Z=12 mm) 



Figure 9. Measured and DSMC Pitot Pressure Profiles for Configuration 1 

(Axial Location. Z = 24 mm) 
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Figure 10. Measured and DSMC Pitot Pressure Profiles for Configuration 1 

(Axial Location, Z = 36 mm) 



Figure 11. Measured and DSMC Pitot Pressure Profiles for Configuration 2 

(Axial Location. Z=0 mm) 
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Figure 14. Measured and DSMC Pitot Pressure Profiles for Configuration 2 

(Axial Location, Z = 36 mm) 



Figure 15. Measured and DSMC Flow Angles for Configuration 2 
(Axial Location, Z = 12 mm) 
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